Anomalous superfluid density in quantum 
critical superconductors 

K. Hashimoto * \ Y. Mizukami * , R. Katsumata * , H. Shishido * ^ M. Yamashita * ^ H. Ikeda * , Y. Matsuda * , 
J. A. Schlueter ^, J. D. Fletcher A. Carrington ' , D. Gnida *% D. Kaczorowski " , and T. Shibauchi * 

* Department of Physics, Kyoto University, Kyoto 606-8502, Japan, ^Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439, U.S.A., I'H. 
H. Wills Physics Laboratory, University of Bristol, Tyndall Avenue, Bristol, BS8 1TL, United Kingdom,* 'Institute of Low Temperature and Structure Research, 
Polish Academy of Sciences, P. O. Box 1410, 50-950 Wroclaw, Poland, * Present address: Institute for Materials Research, Tohoku University, Sendai 980-8577, 
Japan,* Present address: Department of Physics and Electronics, Osaka Prefecture University, Sakai, Osaka 599-8531, Japan, § Present address: RIKEN, Wako, 
Saitama 351-0198, Japan, and ttjo whom correspondence should be addressed. Email: shibauchi@scphys.kyoto-u.ac.jp 

Submitted to Proceedings of the National Academy of Sciences of the United States of America 



When a second-order magnetic ptiase transition is tuned to zero 
temperature by a non-thermal parameter, quantum fluctuations 
are critically enhanced, often leading to the emergence of un- 
conventional superconductivity. In these 'quantum critical' su- 
perconductors it has been widely reported that the normal-state 
properties above the superconducting transition temperature Tc 
often exhibit anomalous non-Fermi liquid behaviors and enhanced 
electron correlations. However, the effect of these strong critical 
fluctuations on the superconducting condensate below Tc is less 
well established. Here we report measurements of the magnetic 
penetration depth in heavy-fermion, iron-pnictide, and organic su- 
perconductors located close to antiferromagnetic quantum critical 
points showing that the superfluid density in these nodal super- 
conductors universally exhibit, unlike the expected T-linear de- 
pendence, an anomalous 3/2 power-law temperature dependence 
over a wide temperature range. We propose that this non-integer 
power-law can be explained if a strong renormalization of effec- 
tive Fermi velocity due to quantum fluctuations occurs only for 
momenta k close to the nodes in the superconducting energy gap 
A{k). We suggest that such 'nodal criticality' may have an impact 
on low-energy properties of quantum critical superconductors. 

superfluid stiffness | d-wave superconductivity | spin fluctuations | mass 
enhancement | quasiparticle scattering 

The physics of materials located close to a quantum critical point 
(QCP) is an important issue because the critical fluctuations as- 
sociated with this point may produce unconventional high tempera- 
ture superconductivity |1. 2|. Quantum oscillations |j3]|4] and spe- 
cific heat measurements |5 | have shown that, in some systems, as the 
material is tuned towards the QCP by controlling an external param- 
eter such as doping, pressure, or magnetic field, the effective mass 
strongly increases due to enhanced correlation effects. Along with 
this the temperature dependence of the resistivity shows a strong de- 
viation from the standard AT^ dependence in the Fermi liquid (FL) 
theory of metals, and often shows an anomalous T-linear behavior 
which corresponds to the A coefficient diverging as zero temperature 
is approached. 

Although there are many studies of non-FL behavior in the nor- 
mal metallic state 1 1 , 2 1, relatively little is known about how the QCP 
affects the superconducting properties below the critical temperature 
Tc. The superconducting dome often develops around the putative 
QCP so that when the temperature is lowered below T^, the super- 
conducting order parameter starts to develop and the Fermi surface 
becomes gapped. It is therefore natural to consider that the low- 
energy quantum critical fluctuations are quenched by the formation of 
the superconducting gap A, which means that the system avoids the 
anomalous singularities associated with the QCP. Perhaps because of 
this reasoning the superconducting properties are usually analyzed by 
the conventional theory without including temperature/field depen- 
dent renormalization effects resulting from the proximity to the QCP. 
For example, in Refs. |6 7| the NMR relaxation rate l/Ti in the su- 
perconducting state is fitted to the temperature dependence expected 



for particular gap functions with the assumption that the normal-state 
1/TiT is virtually temperature independent below Tc even when it 
has strong temperature dependence above Tc due to the magnetic 
fluctuations. 

In superconductors near the QCP, the electron pairing is of- 
ten unconventional with a superconducting energy gap A(fc) which 
changes sign on different parts of the Fermi surface ||8] |9]- This 
sign change stems from a repulsive pairing interaction for example 
resulting from antiferromagnetic spin fluctuations. In many cases, 
this leads to the presence of nodes in the gap A(fe) where the gap 
changes sign. The low-energy excitations from the ground state in 
these superconductors are governed by these nodal regions. In such 
nodal superconductors, the effect of quantum critical fluctuations on 
the excited quasiparticles should be k dependent. Then the question 
arises as to how this effect modifies the low-energy properties in the 
superconducting state. 

The penetration depth A(r) is a fundamental property of the su- 
perconducting state which parameterizes the ability of a supercon- 
ductor to screen an applied field by the diamagnetic response of the 
superconducting condensate. As fermionic quasiparticles are ther- 
mally excited from the condensate a paramagnetic cuiTent is created 
which reduces the screening and increases A. So measurements of 
A(r) give direct information about density and Fermi velocity of 
these quasiparticles |I0|. When the effective mass is enhanced by 
the quantum critical fluctuations, the effective Fermi velocity is ex- 
pected to be suppressed accordingly. In a one-component Galilean 
invariant superfluid, electron correlation effects may not cause the 
renormalization in the low-temperature penetration depth |11|. In 
superconducting materials, however, strong electron correlations do 
affect the renormalization resulting in an enhanced penetration depth, 
which has been reported both theoretically |12 13] and experimen- 
tally 1 141 ITsi . To discuss the energy-dependent effect of quantum 
criticality on superconducting quasiparticles, the temperature depen- 
dence of penetration depth at low temperatures is thus of particular 
importance. 
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Results 

Penetration Depth. Here we begin by presenting results for the 
iieavy-fermion system Ce„nna,i+2 which is located close to a QCP 
(T is the transition metal element, n is the number of Celns layers 
alternating with the Tln2 blocks along the c axis). The most studied 
member of this series is the n — 1 member CeCoIns (T^ = 2.3 K) 
IT6I . in which clear evidence for non-FL behaviors in the normal 
state I17l|18 !| and nodal superconductivity has been found ||6lll9ll20l . 
The recent discovery of superconductivity at ambient pressure in 
the n = 2 member Ce2PdIn8 (Tc = 0.68 K) (21], which exhibits 
very similar non-FL properties [22, 23] |M] |25] [26, 27] to those in 
CeCoIns, allows a detailed comparison of the superconducting prop- 
erties to extract common features in these superconductors near the 
antiferromagnetic QCP. The magnetic penetration depth measure- 
ments have been made for high-quality single crystals of Ce2PdIn8 
and CeCoIns down to ~ 60 mK by using a tunnel diode oscillator 
with the resonant frequency of ~ 13 MHz (see 5/ Text: SI 1). The 
weak ac field is applied along the c axis, which excites supercurrents 
in the ab plane. 

The temperature dependent change in the in-plane penetration 
depth AA(r) = X{T) - A(T = 0) in both CeaPdIng and CeCoIng 
(Fig. lA) exhibits strong temperature variation at low-temperatures, 
much steeper than the flat exponential dependence expected for a 
fully gapped superconductor. The AA(T) data for CeaPdIng are re- 
producible in different crystals and the data for CeCoIns are fully 
consistent with the previous studies |28. 29, 30 1. The strong temper- 
ature dependence indicates substantial excitations of quasiparticles 
at low energies, evidencing the presence of line nodes in the energy 
gap. This is consistent with the residual density of states (DOS) in 
the low temperature limit observed by thermal conductivity 1 19 23), 
specific heat 1201 1241 . and NMR measurements |6 27 |. In particular, 
a d^2_y2 order parameter with nodes along the (110) directions has 
been established in CeCoIns by angle-resolved thermal conductiv- 
ity 1 19 1 and specific heat measurements [20] . The strong similarity 
between CeaPdIng and CeCoIns found in the low-temperature A(r) 
points to common nodal structure in these superconductors. 

In a pure d-wave superconductor with line nodes such as high- 
Tc cuprates, it is well-established that AA(T) shows a linear T- 
dependence at low temperatures, which stems from the linear energy 
dependence of low-energy DOS of quasiparticles. However, a clear 
deviation from the T-linear dependence is observed in the present 
heavy-fermion superconductors (Fig. lA). The data also strongly de- 
viate from the dependence expected for the dirty limit case 



I3II (Fig. 16). We rather find that the power-law dependence T°' 
with an unusual exponent a — 3/2 can describe the observed 
low-temperature variation in a wide range in both superconductors 
(Fig. IQ. 

A few explanations for the super-linear temperature dependence 
of AA(r) in d-wave superconductors have been put forth, including 
the effect of impurity scattering 131], non-local effect near the nodes 
|32l and phase fluctuations |33 |. Often an interpolation formula de- 
scribing a crossover from T to dependence AA(T) oc T'^/{T + 
T* ) is used to describe the experimental data. In particular, the im- 
purity effect leads to the crossover temperature T* ~ 0.83\/rAo, 
where F is the impurity scattering rate and Ao is the maximum gap. 
This successfully accounts for the systematic change of AA(r) with 
impurity scattering observed in Zn-doped YBaaCusOe.gs (Ref. 1 34 1). 
In the present heavy-fermion case, however, fitting to this crossover 
formula in these relatively clean superconductors yields T* /Tc val- 
ues (~ 0.5 for CeaPdIng and ~ 0.2 for CeCoIns) substantially larger 
than the estimates from these theories (see 5/ Text: SI2 and Ref. 1 30|). 
Therefore, the T'^/^ dependence in a wide T/Tc range commonly ob- 
served in these superconductors with quite different Tc, which is dis- 
tinctly different from the T-linear dependence in e.g. YBaaCusOe.gs 
(Ref. 1 34 1), rather suggests some inherent mechanism related to their 
closeness to the antiferromagnetic QCP. 

Superfluid Density. To discuss the precise temperature evolution of 
quasiparticle excitations, we analyze our data in terms of the super- 
fluid density 1/A^(T). We use the reported values of A(0) (280 nm 
for CeCoIns |30| and 1010 nm for CeaPdlng |35|), from which the 
normalized superfluid density paiT) = A^(0)/A^(T) has been ob- 
tained (Fig. 2A). The factor of ~ 3.5 difference in the slope of 
d\/d{T/Tc) (Fig. lA) is consistent with the difference in A(0) in 
these two superconductors, which is also consistent with the larger 7 
value in CeaPdIng |24|. This results in an almost collapse of the full 
temperature dependence psiT) into a single curve (Fig. 2A, inset), 
and the low-temperature variation shows the (T/Tc)^''^ dependence 
with nearly identical slopes. 

This (T/Tc)^/^ dependence can be also seen in other classes 
of materials which are close to antiferromagnetic order. In 
Fig. 2B we show data for the organic superconductor k-(BEDT- 
TTF)2Cu[N(CN)2]Br which is consistent with but measured to lower 
temperature than that in a previous report [|36 :: and the iron-pnictide 
superconductor BaFe2(Aso,7Po,3)2 1151 . Both these materials are 
known to have line nodes in their superconducting gap 1371 1381 . 




Fig. 1. Temperature dependence of the magnetic penetration deptti in heavy-fermion superconductors near the antiferromagnetic QCP. (A): Low-temperature changes 
in the magnetic penetration depth AA(T) = A(T) — A(T = 0) of single crystals of Ce2Pdlng and CeCoIns. The curves are shifted vertically for clarity and the 
data for CeCoIns are multiplied by 3.5. The inset shows the ac susceptibility over the whole temperature range measured by the frequency shift of the tunnel diode 
oscillator, showing sharp superconducting transitions. The dashed line is an exponential temperature dependence expected for a fully gapped s-wave superconductor. 
(B): The same data plotted against (T/Tc)^. The solid line represents a T'^ dependence. (C): The same data plotted against (T/Tc)'^^^. The solid line represents 
a T^/^ dependence. 
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Fig. 2. Universal T^^^ dependence of superfluid density in unconventional superconductors in the vicinity of the antiferromagnetic order. (A): The normalized 
superfluld density ps as a function of (T/Tc)^^^ at low temperatures for Ce2Pdln8 and CeColns. The lines represent T^/^ dependence. The Inset shows 
the overall temperature dependence up to (T/Tc) = 1. (B): A similar plot for Iron-pnictide superconductor BaFe2(Aso.7Po.3)2 and organic superconductor 
K-(ET)2Cu[N(CN)2]Br. The solid lines are the fits to T^/^ dependence. The low-temperature data for BaFe2(Aso.7Po.3)2 are vertically shifted for clarity 



In the BaFe2(Asi_i:Pa;)2 series, there is clear evidence for the an- 
tiferromagnetic QCP being located at a:: ~ 0.30 |15|. In n- 
(ET)2Cu[N(CN)2]Br, although the proposed phase diagram suggests 
that the boundary between the superconducting and antifeiTomag- 
netic states is a first-order phase transition |39 1, the anomalous criti- 
cal exponent near the Mott critical end point |40 | suggests the pres- 
ence of strong antiferromagnetic quantum fluctuations. The normal- 
ized superfluid density ps(T) in these superconductors shows very 
similar [T /Tc)'^^'^ dependence at low temperatures with slight devi- 
ations at the lowest temperatures below T/Tc ~ 0.05 ((T/Tc)'*^^ < 
0.01). 

These results imply that the (T/Tc)^^'^ dependence of superfluid 
density in a wide T/Tc range is a robust property in unconventional 
superconductors, in which strong antiferromagnetic fluctuations are 
present (for comparisons between iron-pnictides and cuprates, see 5/ 
Text: S13). 



Discussion 

If the quantum fluctuations survive the Fermi-surface gapping, the 
effective mass is expected to diverge in the zero temperature limit. 
In such a case, the strong enhancement of mass leads to a reduction 
of superfluid density when approaching the zero-temperature limit. 
It has been recently found in iron-pnictides that the zero-temperature 
superfluid density 1/ (0) shows a strong reduction at the QCP | IS*], 
indicating the strong quantum critical fluctuations directly affect the 
superconducting condensate. However, the temperature dependence 
of the normalized superfluid density ps (T) actually continues to rise 
with decreasing temperature (Figs. 2A and B, insets). This suggests 
that the effect of temperature dependent quantum fluctuations enters 
as corrections in the temperature dependence of superfluid density at 
low temperatures. The anomalous non-integer power-law tempera- 
ture dependence of Ps{T) universally observed in quantum critical 
superconductors thus calls for further theoretical understanding. 

Nodal Quantum Critlcallty. Below we will propose a possible sce- 
nario that this universal behavior can be naturally explained by in- 
voking a strong momentum dependence of renormalization due to 
the nodal gap structure. Suppose that the gap formation below Tc 
quenches the low-energy quantum critical fluctuations, which pre- 
vents the effective mass enhancement in the superconducting state. 
Then in the nodal superconductors, because of the strong momentum 
dependence of the gap magnitude we may consider that the degree of 



quenching of the quantum fluctuations in the superconducting state 
has a strong momentum dependence as well; the low-energy fluctu- 
ations are expected to be strongest near the nodes where the FeiTni 
surface is not gapped (Fig. 3A). To model this effect, we consider the 
angle dependence of the effective Fermi velocity v*p ((j)) along the un- 
derlying FeiTni surface. We consider a simple cylindrical Fermi sur- 
face and assume that the renormalization in Vp{k)/vp is inversely 
related to the enhancement in the effective mass m* (k) (Fig. 36). 
(Here 1/vp is given by the dynamic effective mass which is different 
from, but closely related to the thermodynamic mass II12I ). The ef- 
fective mass enhancement on approaching the QCP can be described 
by m*^ oc (p — Pqcp)~^, where p is a nonthermal parameter con- 
trolling the distance from the QCP at pqcp- The critical exponent 
P has been estimated experimentally by using magnetic fields as the 
parameter p, and in several materials a value of /3 close to unity has 
been reported |41) |42j. In the present case, we take the magnitude 
of the superconducting gap | A| as the control parameter, because the 
gap magnitude determines the degree of quenching of low-energy 
fluctuations. We thus assume vp{k) oc j A(fc)j'^^^, which allows us 
to calculate the temperature dependence of superfluid density by the 
integral over the Fermi surface S as 1 10, 13] 

where Vp (vp) is the effective velocity in the superconducting 
(normal) state, the subscripts i,j denote the directions of the cur- 
rent and vector potential (we take both along a), Y{k,T) — 
-2 — = dE is Yosida function, and f(E) is the 

Fermi-Dirac function for the quasiparticle energy E. By using the 
d^2_y2 formula of A(fc) — Ao cos(2(^) and l3 — 1, the normal- 
ized superfluid density ps is calculated with Vp{(l>) oc | cos(20)|^''^ 
(Fig. 3C), which results in the y^/^ dependence at low temperatures 
(Fig. 3D). 

As real materials will never be located exactly at the QCP there 
will be a cutoff for the diverging m* near the nodes. This can be mod- 
elled by introducing a minimum value for Vp. This leads to an up- 
ward deviation from the t'-^^^ dependence of ps at very low tempera- 
tures approaching T-linear behavior at sufficiently low temperatures 
(Fig. 3D). This can explain the essential features of the experimental 
observations. We note that the presence of the cutoff can be also ex- 
pected even at the exact QCP, because the dynamical susceptibility in 
the zero temperature limit should diverge only on certain regions of 



Footllne Author 



PNAS I Issue Date | Volume | Issue Number | 3 




Fig. 3. Nodal quantum criticality in unconventional superconductors. (A): The momentum dependent gap A(fc) (whose magnitude is illustrated by thin lines with 
grey shades) opens on the Fermi surface (thick line), and has nodes (red circles) at certain directions. In d^2_j,2-wave superconductors, for example, A(fc) has 
strong in-plane anisotropy Ao cos(20) as a function of azimuthal angle (f. In quasi-two dimensional systems, the Fermi surface is approximated by a cylinder, and 
thus the gap has nodal lines at perpendicular to the planes. At the nodes, the gap is zero and thus the quantum critical fluctuations may be present (red shades) on 
the ungapped Fermi surface. (B): The nodal quantum fluctuations lead to the momentum dependence of the renormalization in l/vp{k) near the nodes (blue lines). 
(C): The angle dependence of the renormalized Fermi velocity Vp{(f>) relative to the unrenormalized one vp along the Fermi surface, assumed for calculations of 
the superfluid density in (D). Near the nodes, we illustrate different cutoff levels, which model finite distances from the QCP or disorder. {D): Calculated normalized 
superfluid density as a function of (T/TcY/'^ with different cutoff levels, which explains the deviation from the T'^/'^ dependence at very low temperatures. The inset 
is the full temperature dependence up to (T/Tc) = 1. 



the Fermi surface (determined by the momentum dependence of the 
spin fluctuations), which in general may be different from the nodal 
points |9|. Of course, disorder would also produce some additional 
changes to the temperature dependence but we have not included this 
in the present model. 

An additional factor could also come from the temperature de- 
pendence of the renormalization of m* . As only the theiTnally ex- 
cited quasiparticles will be renormalized, the angular range near the 
nodes where this occurs (i.e., where A(fc) < ^bT) is quite limited 
at low temperature and so this has a rather minor effect. In our tem- 
perature range of interest, T < 0.2Tc, this region is actually limited 
to a narrow angle range of ~ ±3° near the nodes. If we add a cut- 
off constant of 0.3 covering this angle range (Fig. 3C), we found that 
the T-dependence is affected only in the lowest temperature range 
of [T/Tcf-^ < 0.01, above which the T^^^ dependence of ps still 
holds (Fig. 3Z)). This exercise implies that the inclusion of the tem- 
perature dependence of renormalization will not change ps (T) sig- 
nificantly in the temperature range of interest (for more discussion, 
see 5/ Text: SI4). 

We also note that in iron-pnictides the gap symmetry is most 
likely s-wave f43|, and the model based on the d-wave gap may 
not be applicable. However, the fundamental physics that the low- 
energy excitations are governed by the nodal regions should be es- 



sentially the same. Although the detailed structure of the momentum- 
dependent Fermi velocity (such as the precise value of the critical 
exponent /?) will affect the detailed Ps{T) (see 5/ Text: SIS), it is 
striking that such a simple model can capture the salient feature of 
the unusual behavior of quasiparticle excitation in the superconduct- 
ing state of these materials. In a FL theory, the renormalization of 
the effective FeiTni velocity vp in Eq. [1] can be described by the in- 
terplay between fc-dependent quasiparticle interaction and nodal gap 
structure 1 13 1. More detailed theoretical calculations will be needed 
to account for: strong energy dependence of critical fluctuations, the 
effects of disorder (such as quasiparticle scattering interference |44 |) 
and the inelastic quasiparticle scattering rate. 

We suggest that the nodal quantum criticality proposed here is 
an important aspect of unconventional superconductivity close to the 
magnetic QCP. Besides the penetration depth analyzed here there will 
be implications for most other superconducting properties such as 
thermal conductivity and the NMR relaxation rate l/Ti{T) which 
has been long known to have strong deviations from the T"^ law, 
which in the usual analysis gives overestimates of the residual quasi- 
particle DOS in quantum critical superconductors ||6]2]. It should be 
straightforward to extend our analysis to these other properties. 
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SI1-Sample information and experimental techniques 

High-quality single crystals of Ce2PdIn8 and CeCoIns were grown 
by the self-flux method IS1IIS2I . The typical lateral size of the crys- 
tals is 150 X 150 ^ra^ . For Ce2PdIn8 very thin crystals with a thick- 
ness of less than 10 were prepared to avoid contamination of 
Celns. We have checked that our crystals used in this study show 
no antiferromagnetic signals. The magnetic penetration depth mea- 
surements down to very low temperatures (~ 60 mK) have been per- 
formed by using a tunnel diode oscillator with the resonant frequency 
of 13 MHz [S3jlS4], which is mounted on a '^He-'^He dilution refrig- 
erator. The sample is placed on a sapphire rod, the other end of which 
is glued to a copper block on which a Ru02 thermometer is mounted. 
The sample is placed at the center of a solenoid which forms part of 
the resonant tank circuit. A weak ac field is applied along the c-axis 
so that the supercurrent flows within the afe-plane. Changes in the 
resonant frequency are directly proportional to changes in the mag- 
netic penetration depth, AA(T) = GAf{T). The calibration factor 
G is determined from the geometry of the sample |S5 1. 

Samples of K-(BEDT-TTF)2Cu[N(CN)2]Br were grown via the 
standard electrocrystalisation method |S6I. Measurements were also 
made using a tunnel diode oscillator. Special care was taken to cool 
the samples slowly across the temperature region 70 - 90 K where a 
structural ordering transition occurs. The cooling rate was restricted 
to 0.02 K/min, which is sufficiently slow enough to avoid suppression 
in Tc and reduction in electronic contribution to the heat capacity 

(S3. 

The data in Fig. 26 was taken for a sample with dimensions 
0.8 X 0.7 X 0.4 mm^. In these organic samples even weak radio 
frequency (RF) fields can produce self-heating effects at the lowest 
temperatures (T < 300 mK). To eliminate this effect, measurements 
were made with the sample at different positions on the axis of the 
RF coil, effectively reducing the applied field. Typically, a reduction 
in field of a factor three was sufficient to eliminate all signs of self- 
heating. This procedure was tested in three samples and consistent 
results were obtained. 



SI2-Comparisons with existing theories 

The origins of the deviation from the T-linear temperature depen- 
dence of A(r) in a d-wave superconductor are discussed in the frame- 
work of two different theories. Hirschfeld and Goldenfeld |S8 | have 
shown that a small amount of impurity causes a residual density of 
states near Ep, which can change the linear-T dependence (pure 
regime) to a quadratic T-dependence (impurity-dominated regime) 
at a crossover temperature T*^^. This characteristic temperature 
is given by T-*-^p ~ 0.83\/rAo, where F is the impurity scatter- 
ing rate and Ao is the maximum gap. On the other hand, Kosztin 
and Leggett | S9 1 have pointed out that since the coherence length ^ 
diverges near the nodes in a d-wave superconductor, Cooper pairs 
near the nodal directions, where ^ 3> A, cannot participate in the 
field screening. Therefore, the effective penetration depth should be 
larger than the local limit value below a characteristic temperature 
T^i^ ~ Ao^(0)/A(0), which leads to a quadratic T-dependence. 
Thus, for both theories the interpolation formula between two regions 
is given by AA(T) oc T"^ /{T + T*). Here we note an important dis- 
tinction between the temperature dependence of AA(T), which is 
directly measured, and the normalised superfluid density, Ps(T) = 
A^(0)/A^(r). The d-wave form for the superfluid density, ps{T) = 
1 - ^lT/T,, gives AA(T)/A(0) = \KT/T,) + Ip^'iT/T^f + ■ ■ ■, 



1.00 


— «===J 1 1 


1 


0.98 




- 








0.94 






0.92 


• CeCoIng 

• CesPdIng 




0.90 






0.88 


1 1 1 


1 



0.05 0.10 0.15 0.20 0.25 

T/T, 

Fig. S1. Expanded view of normalized superfluid density plotted against T/Tc 
for Ce2Pdlng and CeColns. The solid lines are fit to the crossover behavior in a 
d-wave superconductor described in the text. 

which lead to a small quadratic component to A(r) that depends on 
^, where /j, is proportional to the inverse of the slope of the gap A((^) 
at the nodes, fj, = (4 In 2)fcsrc/(dA(0)/d0|nodc) 1S10|. Then, 
even in the pure d-wave model we have slightly concave tempera- 
ture dependence for A(T). Therefore, we take the following formula 
to describe the crossover behavior expected in these theories; 

p.(T) = i-M(rVTe)/(r + T*). [s-i] 

Figure SI shows the low-temperature superfluid density Ps(T) 
both for Ce2PdIn8 and CeCoIns. We carefully examine the possibil- 
ity of the two theories mentioned above by evaluating T*^^ from the 
present data. The obtained values of T^xp/Tc from the fits are 0.48 
and 0.22 for Ce2PdIn8 and CeCoIns, respectively. If we assume that 
the origin of nonlinearity of A(r) comes from the residual density of 
states induced by impurities, we can estimate the impurity scattering 
rate F via T* ^ 0.83\/rA^. The obtained value of F for Ce2PdIn8 
is O.I 1 K. By using this and Abrikosov-Gor'kov theory|Sl 1 1, we es- 
timated that the Tc reduction is 11% for Ce2PdIn8. However, recent 
specific heat measurements have suggested that the impurity level of 
Ce2PdIn8 is comparable to 0.22% Cd-doped CeCoIns 1S12|. In this 
case, Tc reduction is less than 5 %. Similarly, the obtained value 
of T* for CeCoIns using Eq. [S-1] is much larger than the value 
estimated from thermal transport measurements (< 30 mK) |S13|. 
Therefore, impurity effect are unlikely to explain the nonlinearity of 
A(r) in these systems. As for the nonlocal effects, the reported val- 
ues of ^0 ~ 8nm |S2 S14 S15| for both compounds give T^j, ~ 
14 mK and 140 mK for Ce2PdIn8 and CeCoIns, respectively, which 
indicates that the non-locality effect should be important at much 
lower temperatures than our T^^p- Thus, we can exclude these ef- 
fects as primary origins of the deviation from the T-linear behavior 
mps{T). 
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We also note that the T^/"^ dependence of low-temperature super- 
fluid density has been derived from the Bardeen-Cooper-Schrieffer 
(BCS) to Bose-Einstein condensation (BEC) crossover below Tc 
|S16|. However, this assumes an s-wave state, which is not appli- 
cable to the nodal superconductors focused in the present study. 



SI3-Comparisons between iron-pnicteides and cuprates 

We comment on the difference in the high- Tc cuprate superconduc- 
tors. In clean single crystals of optimally doped YBa2Cu307_4, the 
superfluid density exhibits a rather wide range of T-linear behav- 
ior and the crossover temperature is as low as T*^^/Tc ~ 0.01. In 
hole-doped cuprates, the superconducting dome is separated from 
the antiferromagnetic order, and a possible QCP of the pseudogap 
phase has been discussed near the optimum doping of the supercon- 
ducting dome. The absolute value of the superfluid density shows a 
broad maximum near the optimal doping IIS17I . which is opposite to 
the iron-pnictide BaFe2(Asi_a;Pi:)2 case, where the sharp minimum 
is observed at the antiferromagnetic QCP IS181 . The latter case is 
consistent with the enhanced mass (which enters inversely to the su- 
perfluid density) associated with the QCP. This difference in doping 
dependence seems consistent with the difference in the temperature 
dependence of superfluid density; the y^/^ dependence is a common 
property associated with the mass enhancement due to the antiferro- 
magnetic QCP. 

To see the doping dependence of superfluid density in the iron- 
pnictide BaFe2(Asi_2,Pa;)2, we plot the low-temperature ps as a 
function of {T/TaY'^ in a wide doping range (Fig. S2). The anti- 
ferromagnetic QCP composition in this system has been located at 
X = 0.30 ISI8I . In high composition samples (x > 0.38) which are 
relatively far from the QCP, the low-temperature deviations from the 
T^^^ dependence is more pronounced than those in samples closer 
to the QCP (0.27 < x < 0.33). These results further support that 
the non-integer power law is related to the closeness to the magnetic 
QCP. In hole-doped cuprate superconductors, the antiferromagnetic 
QCP is actually located outside the low-doping end point of the su- 
perconducting dome. It deserves further studies to understand how 
the phase diagrams relate to these differences of low-energy quasi- 
particle excitations in cuprates and iron-based high- Tc superconduc- 
tors. 
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Fig. S2. The normalized superfluid density ps as a function of (T/Tc)^-^ at 
low temperatures for BaFe2(Asi_£Pi:)2 single crystals. The original penetration 
depth data are taken from Ret. IS181 . 

SI4-Temperature dependence of the renormalization 

Near the QCP, the temperature dependence of the renormalization 
is an important factor. In the normal state above the upper critical 
field, the specific heat is largely enhanced at low temperatures |S19 1, 
indicating the strong temperature dependence of m* which may be 
approximated by 1/Vt at low temperatures (Fig. S3A). (Note that 
due to the field-induced QCP, this dependence is the strongest near 
the upper critical field and the actual temperature dependence at zero 
field may be weaker.) In the superconducting state, owing to the mo- 
mentum dependence of the renormalization we consider, l/vp{T) 
will level off at momentum dependent temperature. Near the nodes, 
the temperature dependence of renormalization will continue to low 
temperatures, whereas in the antinodal directions it will become in- 
dependent of temperature near Tc. 

If we assume v*p((j),T)/vF = m&x{C^T/Tc, \ cos(2<j!>)j^''^), 
we can calculate the temperature dependence of averaged mass en- 
hancement in the superconducting state (Fig. S3A). Here we separate 
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Fig. S3. (A): Temperature dependence of the effective mass normalized by its 0.2Tc value In CeColns (black circles) determined from the specific heat 
measured at 5 T in the normal state (above the upper critical field) S19 can be fitted to a l/VT dependence (black line). The inverse of the effective 
Fermi velocity averaged over the Fermi surface In the superconducting state (corresponding to the first term in Eq. [1] In the main text) calculated assuming 
Vp{4>, T)/vp = max(Ci/T/Tc, | cos(2</))|-'/^) with a proper constant C = Vb\ cos 2(45° it 3°)|-'-''^, has a much weaker temperature dependence (red 
line). (B): Temperature dependence of the normalized superfluid density with (solid lines) and without (dashed lines) assuming the temperature dependence of 
renormalization. The total superfluid density (black) can be separated to the diamagnetic (red) and paramagnetic (blue) components. 
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the two terms in Eq. [1] in the main text as A = ^ — Ap ^, 
where the first term A^^ = |^|^ / "i^^f '^^ corresponds to the dia- 
magnetic current and the second term comes from the paramagnetic 
current carried by the thermally excited quasiparticles. In Fig. S3S 
the temperature dependence of each contribution is presented with 
and without the temperature dependence of renormalization. These 
results indicate that in the superconducting state the temperature de- 
pendence of the averaged renormalization is weak, which seems to be 
consistent with the recent high-frequency conductivity results IS20I . 
Moreover, the temperature dependence of the total superfluid density 
A^^ (T) remains almost unchanged even we consider the temperature 
dependence of the renormalization in the superconducting state. 

SI5-Effect of quantum critical exponent 

To see how the power-law of superfluid density varies with the quan- 
tum critical exponent /3, we have computed pa (T) with several values 
of /3 up to ~ 2. The low-temperature 1 - Ps{T) is fitted to the T" 
dependence up to T = Q.2Tc and the obtained power a is plotted 
against the exponent /3 in Fig. S4. We find that the power a varies 
linearly with /3, which can be approximated by a ~ 1 + /3/2. 

The fact that we universally obtained a close to 3/2 in the present 
quantum critical superconductors suggests that the exponent /? in 
these materials is not far from unity. Indeed, in /3-YbAlB4 in which 
zero-field quantum criticality has been suggested IIS2II . thermody- 
namic measurements shows /? ~ 1, which leads to a « 1.54 
(Fig. S4). In the study of YbRh2Si2, where the quantum critical field 
is 0.66 T, /? ~ 1 has also been reported |S22|. We note that in a sim- 
ilar system YbRh2(Sio.95Geo.o5)2 where the critical field is reduced 
to 0.027 T, smaller value of /3 ~ 2/3 has been suggested |S23|, but 
this value gives a ~ 1.36, which is in practise difficult to distinguish 
from the 3/2 power-law dependence. We also note that in CeCoIns 



and Ce2PdIn8 superconductors the estimated values for (3 from trans- 
port coefficient A near the upper critical field Hc2 are rather scattered 
(~ 1.37 for CeCoIns tS24i| and ~ 0.57 for Ce2PdIn8 | S14|), but the 
first-order transition of the low-temperature Hc2 due to the strong 
Pauli paramagnetism P S25I IS 121 may provide some complexity in 
evaluating the precise value of critical exponent. 
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Fig. S4. Low-temperature power a of the superfluid density obtained from ttie 
fitting up to T/Tc = 0.2 for the calculated ps{T/Tc) with different values of 
the quantum critical exponent /3. The line is a fit to a linear dependence. 
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